arXiv:1506.04602vl [cond-mat.other] 15 Jun 2015 


On the Kohn-Sham Approach to Time-Dependent Problems in a 

Density-Functional Framework 

Walter Tarantino* 

Laboratoire des Solides Irradies, Ecole Polytechnique, 

CNRS, CEA/DSM, 91128 Palaiseau, France. 

(Dated: May 2015) 

Abstract 

Predictivity of the Kohn-Sham approach to dynamical problems, when regarded as an initial 
value problem in a time-dependent density functional framework, is analysed for a class of models 
for which the argument devised in the work of Maitra et al. [1] for the standard electronic many- 
body problem does not apply. The original argument is here extended and revised. As a result, 
predictivity for this class of problems seems possible only at the price of introducing extra unknown 
functionals in the corresponding Kohn-Sham equation. Furthermore, the same argument, when 
applied to original electronic problem, suggests that the Hartree-exchange-correlation potential is 
not unambiguously identified by the contemporary and past densities and initial states, but also 
requires knowledge of the divergence of the contemporary Kohn-Sham current. 
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INTRODUCTION 


Time-dependent density functional theory (TDDFT) in the Kohn-Sham approach allows, 
in principle, to calculate the charge density of the electronic many-body problem by solving 
the equation of motion of a system of non-interacting electrons [2, 3]. The same approach has 
been adopted for other many-body Hamiltonians and expectation value of local operators, 
usually refer to as ‘densities’. Predictivity of the original approach was questioned in [4] and 
immediately clarified in [1], 

This work is intended as an extension of the argument of [1] and moves a step forth 
towards a complete description of the mathematical foundations of the Kohn-Sham approach 
within TDDFT as an initial-value problem. 

In the following I will improve the argument of [1] for the standard TDDFT case and 
explain how it does not hold for some other class of many-body problems, including some 
models of physical interest. I will discuss the consequences of the failure of such an argument, 
showing how the issue of ‘predictivity’ can be discussed in terms of existence of a unique 
solution of the equations of motion of the Kohn-Sham system. I shall then characterize such 
a class of problems and, finally, explain how in some cases the approach can still be made 
‘predictive’, as long as new (unknown) density-functionals are introduced in the relevant 
equations. 

In Section I the Kohn-Sham approach within the framework of Density-Functional Theory 
for time-dependent problems is recalled; in Section II a simple model (Hubbard dimer) for 
which such a procedure is problematic is presented; in Section III a tool for understanding 
whether a model suffers of the same problem is devised and a simple modification of the 
original recipe that fixes the problem for the Hubbard case is generalized. 


I. THE DENSITY-FUNCTIONAL TIME-DEPENDENT KOHN-SHAM APPROACH! 
(1 OF 2) 

TDDFT in its KS declination represents an alternative to the time-dependent Schrodinger 
equation for describing the evolution of the charge density of a system of interacting electrons 
[5]. In fact such a reformulation of the electronic problem can be generalized and used for 
other time-dependent problems as well. The common way to perform such a generalization 
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and apply it to other models is the object of this section. 

We start with a generic time-evolution quantum problem, formulated in second quanti¬ 
zation and characterized by the equation of motion 

id t \y(t))=(f + W + V(t^\y(t)), (1) 

which includes a kinetic term T, an interaction term W and the action of an external classical 
held via V(t), and by a given state |'h 0 ) such that 

|*(« = 0 )) = |* 0 ). ( 2 ) 

I assume that the first order differential equation (1) has one and only one solution | *&(£)) 
satisfying (2). I shall refer to this as the ‘time-dependent many-body problem’, or, simply, 
‘the many-body (MB) problem’. 

If we are not interested in finding the state vector that solves the problem, but only in 
the expectation values of a finite number of local operators (from now on ‘densities’), we 
can attempt to approach the problem via the Kohn-Sham method, which consists in solving 
an auxiliary differential equation that is only linear in the held and characterized by a set 
of effective potentials such that the densities evaluated in such a system match their value 
in the original problem. Let us write V(t) = XWext (t)n, where the sum is a shorthand 
notation for a integration/sum over all residual degrees of freedom (space/lattice points, 
vectorial/tensorial/spin components...) and n is a local operator bilinear in the held. For 
sake of argument, I consider the case in which we are interested only in the expectation 
value of n, namely {('F(t)|h| v L(t)) | t 6 [0, +oo)}. The Konh-Sham system for this problem 
is then identified by the equation of motion 

a, |tw> = (f +k«) i»«> (3) 

with V s (t) = J^ y s (t)n, and a certain initial state 

|$(t = 0)) = |$ 0 ) (4) 

such that ('F(t)|h|'F(t)) = ( < f>(t)|h| < h(t)) for all t e [0, +oo). The existence of such a Kohn- 
Sham system is in general not known and it becomes crucial to prove it even under restrictive 
assumptions. To proceed with my argument, I assume that a theorem similar to the van 
Leeuwen-Runge-Gross theorem [2, 3] can in fact be proven for the model under exam. Such 
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a theorem ensures existence and uniqueness of the system, under specific conditions, and 
allows to rewrite the problem in a density-functional framework. It can be stated in the 
following terms. 

Given f(t) a generic function of time, we distinguish the set {/(f) | t e [0, +oo)} de¬ 
noted by simply /, from the value of the function at the specific time t denoted by f(t). 
A set of functions is then denoted by $ = {f,g, ...}. The expectation value of an observ¬ 
able O is denoted by 0 mb( 4 if sandwiched by |T(t)), or by Oics(t), if sandwiched by 
|<L(f)). Moreover, n(t ) = umb(^) = riKsit). Then, given a certain set of external potentials 
•(Text) v' xt , V "xti we denote by 91 = {n, n', n", ...} the set of densities obtained by 

solving the problem (1-2) with any element of 9J C (giving rise to the set of ‘v-representable 
densities’). The first part of the theorem is then stated as following: 

Theorem 1.1 There is a one-to-one correspondence between elements o/9J c and elements 
of 91, up to gauge transformations. 

Moreover, given the set 911 of densities obtained by solving (3-4) by means of some effective 
potentials, the second part of the theorem reads 

Theorem 1.2 There exists a set 9/ = {u s , v' s , v", ...} such that 931 = 91 and, moreover, 
elements of%3 s are in one-to-one correspondence with 91, up to gauge transformations. 

From now on I shall refer to these two statements as to ‘vLRGl’ (van Leeuwen-Runge-Gross- 
like theorem). For the electronic many-body problem the theorem has been proved under 
some restrictive assumptions [2, 3, 6], which for the moment are not of our concern. 

Once the existence of the Kohn-Sham system is ensured by vLRGl, one has to find a 
viable way to calculate the effective potential v s . Looking back at the inputs of the original 
problem, v s (t) can be regarded as a functional of |T 0 ), |$o) an v ext - This will be denoted 
by v s (t) = F t [u ext ], in which the explicit dependence on initial states is omitted but always 
understood. If we were given this functional, then the problem (3-4) would have been well 
defined and admit, like the problem (1-2), one and only one solution. 

A density-functional formulation of the problem, which we are allowed to do thanks 
to vLRGl, implies the definition of quantities that depend on the density and not on the 
external potential. We cannot use directly v s (t) = F t [n] , even though vLRGl ensures its 
existence and uniqueness, since the problem (3-4) would have no information at all about 
the external potential that drives the evolution of (1). 
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To feed in some information about the external potential, one usually [2] writes 

V s (t) = V ext (t) + UHxc (t) (5) 

which in fact defines a new quantity, vu xc (t). As long as vLRGl holds, we can write uhxc(^) = 
G t [n], which expresses the functional dependence of vn xc (t) on n. 

It should be emphasized that vLRGl states that, given the entire set {n(£)| t e [0, +oo)}, 
the entire set {t'Hxc(^)| t G [0,+oo)} is unambiguously identified; however, it does not 
specify how single elements of these sets are linked. More specifically, it does not tell 
whether information on future densities is necessary to construct t’Hxc(f) or not. As discussed 
in [1, 4], it appears that this might threaten the predictivity of the approach. For the 
standard electronic many-body problem of [2], an argument to prove that these worries are 
unjustified was presented in [1], in which it was explained how r?Hxc(£) requires no information 
on the future to be completely determined. However, this argument is specific for the 
Hamiltonian and the choice of densities of [2], Moreover, part of it is in conflict with the 
initial assumptions of [2], as I shall soon explain. In the following section I shall present a 
simpler model, which will allows me to enlighten some limits of the original argument of [1] 
and suggest possible solutions. 


II. THE CASE OF THE HUBBARD DIMER 

A. Definition of the model and corresponding Kohn-Sham system 

I now consider the Hubbard dimer in presence of an external field. Such a model has 
been analyzed in a TDDFT framework in several works [7-9] but its relevance to the realistic 
electronic many-body problem is not of our concern, as I use it only as prototype of a class 
of problems for which the argument of [1] does not apply. The model is defined by 

H = (cLc/fe + cLdLa) + U (' h L fh Li + n^h^) + v(t) — ^ (6) 

<T=t,t 

where h t = XLfLuct, and {c lcr , cj p } = 8ijd ap , with i,j = R, L and er, p =j\|, all other 
anticommutators being zero, and by a given initial state |T(f = 0)) = |To). 

If one is interested only in the value of (fix,) and (h r), one can attempt to solve the 
problem by means of the approach above outlined. In fact, the value of (fix, + hu) is only 
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determined by the initial state, since [hi + ha, H(t )] = 0. We can then concentrate on the 
observable n(t ), where 

h L — h R 

n = —2~- (7) 
Following [9], I chose the corresponding Kohn-Sham to be characterized by the following 
Hamiltonian: 

H s = -rY^ {c[ a c Ra + + v s (t)h. (8) 

o-=t,F 

For simplicity, I consider |$(f = 0)) = | H/o) - Then, following the procedure outlined in the 
previous section, I define 

^Hxc (t) =v s (t)-v(t). (9) 

Such a function is determined by the condition 

= riKsit)- ( 10 ) 

This condition can be recast as 

' 

n-Ksi 0 ) = % b ( 0 ) 

< h,Ks(0) = hMs(0) (11) 

nKs(t) = n M B(t), 

The hrst two identities are satisfied no matter what effective potential is chosen, as long as 
the two initial states are taken to be the same. So we can say that vn xc (t) is in fact defined 
by the condition 

nRs(t) = h M B(t). ( 12 ) 

The quantity h(t) can be connected with the contemporary external potential by means of 
the equations of motion (3) and (1): 

^ext(t) (T) MB = h(t) + 4r 2 n(f) + (p, T\,W]) MB . (13) 


Similarly, for the KS system we have: 

v a (t) (T) KS = h(t) + 4r 2 n(f). 
In fact, using (12,13,14,9) it is possible to derive that 


VHxc (t) = V(t) 


(T) 


MB 


([[h,T],W]) 


MB 


(T) 


KS 


(T) 


KS 


Vice versa, as long as (15) holds, equation (12) and hence (10) follow. 


(14) 


( 15 ) 
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B. Taylor-expandability and Causality 


Adapting the constructive algorithm of [3], it is possible to prove that a set of uh X c 
satisfying (15) does exist, at least if we restrict our attention to the set of external potentials 
that are Taylor-expandable in the time variable (‘Te’ potentials, from now on). In fact, under 
this restriction, adapting the entire argument of [3] allows to prove the vLRGl for this case. 
This also allows us to density-functionalize the theory. We can then correctly regard t’Hxc(t) 
as a functional of the density only. I now want to investigate the property of ‘predictivity’ 
of the approach, in the sense of uniqueness of the solution of the KS system (3,4), with 
vhxc (t) = G t [n\. 


Following [1], one might start by wondering about the dependence of t'Hxc(^) on preceding 
(0 < t' < t), contemporary (f = t), and successive (t/ > t) densities n(t'). In [1] it was 
argued that, for any T > 0, vn xc (t) with 0 < t < T cannot depend on densities n(t) with 
t > T, using the Runge-Gross theorem and the fact that two external potentials that are 
the same in [0, T] but differ for t > T must give different densities at times t > T. However, 
for Te external potentials this very first step is problematic, for two Te functions equal 
in a finite interval 0 <ei<t<e 2 <T are in fact equal everywhere within the radius 
of convergence of the Taylor series. This means that, when Taylor expandable external 
potentials are considered, as we did in order to prove the vLRGl, two external potentials 
that were the same in a finite interval in the past will be necessarily the same also in the 
future (up to the radius of convergence, i.e. until vLRGl holds). Equivalently, within the 
radius of convergence of the Taylor expansion of the density, any dependence on the density 
on a finite interval in the past can be recast as a dependence on the density on a finite 
interval in the future, and vice versa. In this case, the intuitive notion of t'Hxc(t) being 
‘causal’, in the sense that does not depend on future densities, becomes meaningless. 


Nonetheless, a meaningful concept of ‘causality’ is not necessary to discuss the predictivity 
of the approach, that remains a non-trivial problem. In order to proceed, we look back at 
equation (15) and consider separately the objects in the form of (O) from v(t). 
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C. Density Functionalization (lof2) 


The contemporary expectation values (O) in (15) are unambiguously identified by the 
state vectors |\I/(£)) and |$(t)). In order to argue that these objects, when regarded as 
functionals of the density only, pose no threat to the time propagation of the Kohn-Sham 
system, we extend the second part of the argument of [1], based on a propagation on a finite 
time-grid. 

More specifically, we consider a discrete time variable and a simple rule for translating a 
differentiation on a continuous variable, namely 

d f(t) fi+i-fi ' fl6) 

dt A 

To compare results on the time-grid with the continuous limit, we consider only terms of 
leading order in A, since the limit A —» 0 makes (16) an identity. 

If we denote the many-body and Kohn-Sham state at the time step i with %p t and 0*, 
respectively, we can write 


A+i = F(vi,ipi) (17) 

<Pi+i = G(vf , (f>i) (18) 

Vi = V(n i+ 2 , n i+ 1 , ipi, (pi) (19) 

vt = V s MiAi) ( 20 ) 

n i+ 1 = N(n h Ji) (21) 


with F(.) obtained from (1,6,16), G(.) from (3,8,16), V(.) from (13A6), 1 V s (.) from (9,15), 
and N(.) from the continuity equation 

n(t) = -m)\J\m), ( 22 ) 


with J = i[h,T]. Then, for the two terms appearing in (15) we can write 


for which 


(Gi) mb 

m KS 


->■ (pi) 


(23) 


1 In fact, an entry in V(.) for (pi is redundant, since Vi is only determined by n,:+ 2 , rij+i and %pi. However, 
this is specific for the model here considered, while in general also (pi may be required. 



L{Mi) = L(F(vi-i, G(v ? s _i, ipi-i)) = (24) 

= L(F(V (rij+i, rii, V’i-i, <f>i-i), V’i-i) ) G (V s (V(ni + i,rii, V’i-i, 0i-i), V’t-i) 0»-i)) V’i-i)) = (25) 

= M(n i+ i, n*, V’i-i, 0*-i) = (26) 

= M(N(ni, Ji),rii, tpi-i, <j>i~i) = (27) 

= f^ii ifti—li 0i—l) (28) 

= I/ ( i) (Jj, rii, F(vi- 2 , tpi- 2 ), G{yl_ 2 , ipi- 2 )) = (29) 

= L (1) ( Jj, ni, F (V (n i} rii _i, ipi- 2, pi- 2 ) , ^-2), <3 (K S (K(rij, 7 ij_i, ^-2, Pi- 2), ^-2, Pi- 2 ),Pi- 2 \Wt 

= -7/(2) ( Ji? Tli—li Pi—2 1 (pi— 2 ) (31) 

; (32) 

■7'(®) ( j 1 ? • • • 7 ^ 17 'Po 7 00 ) • (33) 


This means that on a time-grid these terms are unambiguously determined by: the Kohn- 
Sharn state vector at same time step (which enters through Jj), the set of contemporary 
and previous densities, and the initial states. Since at a given time step, these ingredients 
are in principle all at our disposal, they pose no threat to the time-propagation of such a 
discretized KS system. While rigorous for the discretized problem, this argument is only a 
first step towards a complete proof for the continuous time problem. 

Nonetheless, when applied to the electronic problem of [2], the argument supports the 
claim of [1] that terms in the form of (23) pose no threat to the time-propagation of the 
corresponding KS system. 

In fact, following farther the argument, we are also led to conclude that the rH XC (r, f) 
of [2] cannot be determined by sole knowledge of initial states, contemporary and past 
densities, but also requires knowledge of the divergence of the contemporary Kohn-Sham 
current V • j(f), which appears in the equivalent of (22). The importance of such a term 
was recognized in [1], but only at t — 0, when its value can be obtained from the initial KS 
state. Pushing farther their own argument, as done above, suggests that the dependence of 
UHxc on V • j for finite times, rather than remaining a dependence on the initial V • j (t = 0), 
it becomes a dependence on the contemporary V • j. 

Such a dependence is not incompatible with vLRGl, for which we expect that complete 
knowledge of the density suffices to identify the system. Jj has indeed been introduced to 
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encode the information contained in rq+ 1 , which, in the language of the continuous time case, 
translates into hit). Choosing of expressing such an information in terms of J(t), which one 
can express in terms of |<f>(t)), rather than h(t ), avoids the complications of having a time 
derivative of the state on the right-hand side of the equation of motion (3). 

D. Density Functionalization (2of2) 

Once the n dependence of the terms ( O) has been clarified, to some extend at least, we 
can go back to (15) and discuss the presence of v(t). This needed to not to be discussed in 
[1], because of no explicit v(t) dependence in the corresponding equation. More precisely, in 
the standard TDDFT case, the equivalent of (13) and (14) reads 

V • [n(r, t)Vv(r, t)] + (34) 

and 

V • [n(r,t)Vv a (r,t)] + (0 2 ) MB (35) 

with 0 2 and 0\ two appropriate operators, n(r,t) the electronic charge-density, v(r,t) the 
external held and v s (r,t) the effective held of the corresponding KS system [3]. Since both 
potentials are multiplied by the electronic charge-density, a quantity that is the same for the 
many-body and the KS system, the explicit u(f)-dependence in the corresponding v Bxc (r,t) 
cancels out. 

The fact that in general (T) MB / (T) KS does not simplify to 1, leading in our case to 
the missed cancellation, can be proved by considering the specific case of t — U — 1, 
v(t) = sin(t), and |To) being the half-filled ground-state of the Hamiltonian at t — 0, for 
which 

(?)mb = ~ 2 + *-ll*+$Zt* + 0(1») 

(f) KS = -2 + ^ + %G + 0(t*), 

as one can prove by using the constructive algorithm of [3]. 

The expression on the right-hand of (15) side is therefore not ‘universal’, in the sense 
that it explicitly depends on the system via v(t). However, this dependence can be removed 
if we use again equation (13) to extract the information about the density n contained in 
v(t). This leads to 

v(t) = F t [n\ + G t [n\h. (37) 
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with F and G some suitable functionals. While for these two functionals the argument of 
the previous section applies, the presence of hit) needs a dedicated analysis. 

Even though one could discuss the consequences of the presence of hit) on the time- 
grid, in fact we do not need to, as they can be cleared enough by directly looking at the 
original problem in continuous time. When (37) is plugged back into the equation of motion 
(3), an explicit dependence on the second-time derivative of |<h(f)) is introduced. Even 
assuming that the conclusions of the previous sections can be safely generalised to the 
continuous time case, the presence of a double time-derivative of |<h(f)) changes the status 
of the equation (3) from Erst to second-order. This implies that the initial value (4) is no 
longer sufficient to ensure one single solution to the problem. In other words, the density- 
functionalization of the Kohn-Sham system (8) leads to a problem that, despite being an 
actual exact reformulation of the many-body problem, does not admit only one solution, and 
hence it fails to unambiguously characterize the solution of the initial many-body problem, 
even if all functionals were known. Using the terminology of [1, 4] one could say that this 
Kohn-Sham system is not predictive. 


E. Restoring Predictivity 


Applying the recipe of Section I to the Hubbard dimer above defined led to a reformulation 
of the original problem that, despite being exact, cannot be used to calculate the wanted 
densities. 

This problem can in fact be avoided if we use a slightly modified recipe, which however 
has the cost of introducing a new (unknown) functional of the density. More specifically, 
the way information about v ext (f) is extracted from v s (t) has to be modified as follows: 


v a (t) 

< Vi(t) 

V2 (t) 


Vi(t)v ext (t) +v 2 (t) 


F) mb 
{T) ks 

{[[n,T],W]) MB 
(Fks ' 


(38) 


When Vi and v 2 are regarded as functionals of n, the argument used above for (O) applies 
again and the problems cause by the presence of h(t) are avoided. 
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III. THE DENSITY-FUNCTIONAL TIME-DEPENDENT KOHN-SHAM APPROACH! 
(2 OF 2) 

A. A New Recipe 

The case of the Hubbard dimer here considered is not an isolated exception, as we shall 
see in the next section. In fact we can use what learned from this example to identify a class 
of problems that, if approached in the standard way, present the same complications. 

Given an Hamiltonian and a certain choice of densities for which the vLRGl holds, we can 
connect the external potential(s) to the effective one(s) by differentiating the densities and 
using the equations of motion (1-3). This would lead to an expression that in the simplified 
notation of Section I would read 

{Oi [Us(£)])a-s + (0-2)KS = (®3[Vext(t)])MB + (^4 )MB (39) 

with Oi[f(t)} acting as an operator on the Hilbert space of state vectors and possibly a 
differential operator on the argument fit). If 

mfm KS ~ mfm 

MB > (40) 

then we can say that the problem (Hamiltonian+choice of densities) is RG-like and the KS 
system is as predictive as the one of standard TDDFT [2]; otherwise, we can say that it is 
RG-unlike and considering the Hartree-exchange-correlation potential defined by uhxc(^) = 
v s {t) — v ex t(f) as a functional of the density only leads to a reformulation of the many- 
body problem that does not admit one unique solution, making the density-functionalized 
Kohn-Sham system not predictive. 

In case the operator Offit)] is purely multiplicative on its argument: Offit)] —> Oif(t), 
like in our example of the Hubbard dimer, it is possible to fix this problem by simply 
considering 

' 

V s (t) = 'UlW'Uext (t) + V 2 (t) 

< = (41) 

\yi Iks 

l H(t) = mz 

and regarding V\ (t) and v 2 (t) as functionals of the density separately. 

When the density(/ies) is(/are) chosen to be the expectation value of the local operator(s) 
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conjugate to the external potential(s) in the many-body Hamiltonian, one has 


(6, Mt)]> = ([b,[f,«]]) 

(42) 

<0 2 ) = <ff, [f, ft]]) 

(43) 


(44) 

(0 4 > = ([T + W, [f+W,h]]). 

(45) 


by virtue of [V, n] = 0. 


B. Relativistic Electrons 


Another example of the class of RG-unlike problems is provided by the Kohn-Sham 
approach to QED defined in [10], to which 1 refer the reader for the notation of this section. 
In that work the four-potential is chosen as ‘density’. As one can easily verify from 
equation (31) of [10], our equation (39) becomes 

(0^ U ) KS a t + (Qkin)KS = (0 ^“ / )qed c C Xt + (2km + Shit )QED (46) 


with 

6^ = $ - 7*7 V) '•P (47) 

where ^ is the Dirac held operator, t/5 = t/dy 0 and y^ the gamma matrices acting on the 
spin degrees of freedom of and in which explicit spacetime dependence has been omitted. 
Since (0^ v ) KS ^ {0 ^ v )mb the approach has the problem of predictivity here discussed, if 
a Hxc( a; ) i s regarded as a functional of the four-current only. 2 

The case of QED is particularly instructive because it shows that the predictivity issue 
here discussed does not seem to be necessarily inherent the many-body Hamiltonian con¬ 
sidered, for the choice of densities and the Kohn-Sham system can also be crucial. In [11] a 
different set of densities was chosen for the same QED problem, namely the four quantities 
^dy^ = PA In this case, the equivalent of (39) calculated in the Coulomb gauge (using 
2 In applying the argument of Section IID one might think that using again equation (31) of [10] to connect 
a ext( x ) to p, the resulting Dirac-Kohn-Sham equation would not become of second order, but remaining 


of first. If so, the argument on the number of solutions would no longer apply. However, only two 
components of a[f xt ( x) are determined by (31), since the kernel of the matrix (CP") is of dimension 2. 


Even if one component is fixed by a gauge condition, another component remains to fix and for this one 


should necessarily look at higher derivatives. 
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equations (63) and (65) of [11]) leads to (Oi[v s (t)}) KS — > P°a s and (Os[v ex t(t)]) MB —> P°a ext . 
Being P° one of the chosen densities, the identity ((Oifuext^)])/^ = (C^fext (t)}) MB is i n fact 
fulhlled making the system of [11] RG-like. 


CONCLUSIONS 

In this work the arguments used in [1] to support the predictive character of the Kohn- 
Sham system as an initial-value problem within a TDDFT framework were reviewed and 
extended. 

It was pointed out that the original argument cannot be applied if external potentials that 
are Taylor expandable in time are considered. Then, in order to study the ‘predictivity’ of a 
density-functionalized Kohn-Sham system, identified with the property of such a system of 
admitting one and only one solution when regarded as an initial-value problem, an argument 
on a time-grid of infinitesimal pace has been developed. When applied to the standard 
TDDFT case, such an argument, on one hand, confirms the claims of [1] about predictivity 
of standard TDDFT; on the other hand, it suggests that, in order to preserve the differential 
structure of the Kohn-Sham equation and have a predictive Kohn-Sham system, uhxc(u t) 
must be regarded as functional of contemporary and previous densities, the many-body and 
Kohn-Sham initial state and the divergence of the contemporary Kohn-Sham current. The 
argument is however rigorous only on a time-grid and the problem demands for further 
investigations. 

The specificity of the original argument of [1] to the TDDFT problem has also been 
enlightened by showing another many-body problem, a time-dependent Hubbard dimer, 
whose Kohn-Sham system is in fact not predictive, as long as the corresponding t'Hxc(f) = 
v s (t) — v ext (t) is regarded as a functional of the density only. It was proved that the 
Kohn-Sham equation resulting from a density-functionalization of uhxc was requiring more 
boundary-conditions then provided, failing to being characterized by a unique solution, and 
hence failing to make a ‘predictive’ system. The origin of such a problem was identified, 
allowing to derive a criteria to fulfill for a Kohn-Sham system for not being affected by the 
same predictivity issue. 

In some simple cases, like for the Hubbard dimer here considered, predictivity can be 
restored by modifying the definition of the quantities to be considered functionals of the 
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density. More specifically, it was argued that, if the effective potential is decomposed as 
v s (t) = vi(t)v ext (t) + V 2 (t), with vi(t) and v 2 (t) unambiguously defined as expectation values 
of some specific operators, and the two potentials v\ (t) and v 2 (t) are regarded as functionals 
of the density, the Kohn-Sham initial-value problem has in fact only one solution. 
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